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Considering the propagation of fields in the spacetime continuum and the well-defined features of fields with finite 
degrees of freedom, the wave function is expanded in terms of a finite set of basis functions localized in spacetime. 
This paper presents the analytic eigenenergies derived for a confined fundamental fermion-antifermion pair under 
a linear potential obtained from the Wilson loop for the non-Abelian Yang-Mills field. The Hamiltonian matrix 
of the Dirac equation is analytically diagonalized using basis functions localized in spacetime. The squared lowest 
eigenenergy (as a function of the relativistic quantum number when the rotational energy is large compared to the 
composite particle masses) is proportional to the string tension and the absolute value of the Dirac’s relativistic 
quantum number related to the total angular momentum, consistent with the expectation. 


1 Introduction 

In the formalism elaborated in our previous publications 
mm, the fields are expanded in terms of basis functions 
localized in spacetime. A key characteristic of this for¬ 
mulation based on finite element theory is that it is pos¬ 
sible to apply differentiation unlike in the finite difference 
method. Our method allows the use of a basis set of step 
functions, which is rather different from the formulation 
by Bender et al. [3] . In the non-Abelian Yang-Mills case 
[4-7], the analytic continuum classical field as a possible 
vacuum reveals the linear potential, and quantum fluctu¬ 
ations are expressed in terms of step functions exhibiting 
Coulomb potential. Regarding the confined bound state 
of a fundamental fermion-antifermion pair and related 
themes, one can observe approaches from field theory and 
lattice gauge theory in prior literature [8-18]. Although 
the Dirac equation with a linear potential was investi¬ 
gated analytically by other authors [19-22], the Regge 
trajectory [23] may not be reproduced systematically in 
the relativistic scheme. Concerning a classical mechani¬ 
cal Hamiltonian that describes the principal properties of 
the Regge trajectory, comprising a linear potential and 
repulsive rotational potential [23] [23], which is to be de¬ 
scribed in the part of this paper containing Eqs. J57D- 
(l62l) of Subsection 3.1, the Hamiltonian has no basis in 
the theory of Dirac fields with the potential produced by 


the Yang-Mills fields. Furthermore, a clear answer has 
not been fully provided at the quantum level by other 
theoretical/numerical approaches to the following ques¬ 
tions. What is the mechanism of the mass of a pair com¬ 
posed of a fundamental fermion and antifermion? Why 
can the binding energy (mass) of the paired fermions be 
large compared with the masses of the composite fermion 
and antifermion? The confined bound state of composite 
fermions is not fully understood theoretically. Because 
the Regge trajectory is expressed by the quantum number 
corresponding to the rotational quantity, it is expected 
that the Dirac fields, expressed in spherical coordinates 
with the potential derived directly from the Yang-Mills 
theory, will provide answer to the aforementioned ques¬ 
tions and reveal the meaning of the classical mechanical 
Hamiltonian. 

This work is aimed at obtaining an analytic expres¬ 
sion for the eigenenergies of a confined fundamental 
fermion-antifermion pair using our formalism mentioned 
above mm, which formulates fields of finite degrees of 
freedom using basis functions localized in the spacetime 
continuum. The present formalism enables an analyti¬ 
cal calculation without using numerical values unlike nu¬ 
merical computer simulations via numerical values. The 
action-like total Hamiltonian including a given linear po¬ 
tential, which leads to the Dirac equation in spherical 
coordinates by variational calculus, is expressed in terms 


1 


of a basis set of step functions localized in spacetime. 
The Hamiltonian matrix in the secular equation is diago¬ 
nalized analytically, and the lowest eigenvalue is derived 
as a function of the string tension and Dirac’s relativistic 
quantum number related to the total angular momentum. 
The squared system energies for the large rotational en¬ 
ergy compared to the constituent particle masses corre¬ 
spond to those in the classical mechanical Hamiltonian 
case, which are consistent with the principal properties 
of the Regge trajectory [23]. We emphasize that (1) the 
analytical approach shows clearly that the squared sys¬ 
tem energies are proportional to the string tension and 
quantum number, and that they originate from the secu¬ 
lar equation structure leading to the absolute value of the 
quantum number that includes a sign and constructs the 
potential; (2) the quantum number, instead of being a 
non-relativistic value, is the relativistic integer described 
by the Dirac equation for a pair of fundamental fermions. 

This article is organized as follows. Section 2 de¬ 
scribes the formalism for the Dirac equation using a ba¬ 
sis set of step functions localized in spacetime. Section 3 
presents an analytic expression for the eigenenergies of a 
confined bound fermion-antifernrion pair, including dis¬ 
cussions in Subsection 3.2, and is followed by conclusions 
in Section 4. 

2 Basis equations and theoretical 
formalism 


(—m — —)G + — —F = EG, (5) 

r dr r 

where m is the fermion mass and E is the fermion energy. 

Considering the forms of i/>f = F/r and ipc = G/r 
as well as the two-dimensional integral f dr r 2 , the total 
Hamiltonian that variationally leads to the above Dirac 
equations is 

U = \ f dr[F(+m - —)F - F^- - F-G 
2 J r dr r 

ac _dF 

+ G(—to- )G + G— - G-F 

r dr r 

— FEF — GEG]. (6) 

We now add the linear energy potential derived from 
the Wilson loop for a non-Abelian field. Energy is one 
component of the four-vector momentum and the lin¬ 
ear potential constructs the energy potential combined 
consistently with the Coulomb potential, the Coulomb 
potential —a/r above is replaced by —a/r + or. Here, 
a = g 2 /(4n) with a coupling constant g, and a refers to 
the string tension. We then have the total Hamiltonian 
for the non-Abelian case 

H = ^ [ dr[F(+m — — + ar)F — F — F—G 
2 J r dr r 

a „dF 

+ G(-m -b ar)G + G— - G-F 

r dr r 


2.1 Dirac equation for the non-Abelian 
case in spherical coordinates 

To avoid confusion, we first note that the relativistic n 
used in this paper has the following relation [26] with the 
Dirac’s notation jn m 


- FEF -GEG). 


(7) 


K = JD = 


l + l for j = l + 1/2 
—I for j = l — 1/2 


(1) 


where l and j refer to the quantum numbers for angular 
and total angular momentums, respectively. For the low¬ 
est energy case of a hydrogen atom, n > 0. Denoting the 
radial wave functions of fermions in spherical coordinates 
with radial r-axis as 


2.2 Fields expanded in terms of basis 
functions localized in spacetime 

The present formalism for fields is firstly based on the 
propagation of fields in the spacetime continuum, and 
secondly, on the fact that fields are definite in a scheme 
of finite degrees of freedom. We then expand the fields in 
terms of basis functions localized in the spacetime contin¬ 
uum, which has a finite number of lattice (grid) points, 
by realizing the following formulation. For the consider¬ 
ing region in spherical coordinates, we introduce lattice 
(grid) points r n (n = 1,2,...., N) in the radial r-axis. The 
infinitesimal positive A is defined by A = r n+ 1 — r n , and 
r n ~ 1/2 = r„-(l/2)A and r n+1/2 =r n + ( 1/2) A. The ba- 


^r )= F{T \ 

r 

(2) 

sis functions, which have a superscript without and with 
prime, are defined by 

i f ^ G w 

^G{r) = 

r 

(3) 

fi ” (r) = { 0 

for r n _i/ 2 < r < r n+1/2 

. (8) 

for r < r n _i/ 2 or r > r n+1/2 

the Dirac equations in natural units become 

(+to - -)F G = EF , 

r dr r 

(4) 

(r) = | ^ 

for r n _i/ 2 < r < r n+1/2 

, (9) 

for r < r n _i/ 2 or r > r n+1/2 
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which have the derivatives 

dr \r=r n -i/2 ~ dr lr=r n - 1/2 

= S(r-r n _ i/ 2 ), (10) 

«r), _ <'(>•), 

dr k=r " +1/2 _ dr |r=r " +1 / 2 

=-<5(r-r„ +1/2 ). (11) 

The above basis functions have the following proper¬ 
ties 


\ [ dr(n°n f + nf nf + + nf nf) 


^d nk A — A8 n ki 


(18) 


drll^(r) 


dr 


■ /r>B df2jF 

dr{n ‘nt 


n dr 


n dr 


n dr ■ 


— -[0 + 0 + (i5fc jn+ i — dfc jT i_i) + (<5fc jn -)_i — <5fc,n-i)] 


J drfl f (r)fif (r) = J drfl% (r) = Ad nfc , (12) 

J drfl f (r)flf (r) = <5 nfe J drfl% (r) = Ad nfe , (13) 
J drfln (r)nf (r) = S nk J drfl%'(r) = A S nkl (14) 

where <5 ra fc is Kronecker’s delta, and considering the over¬ 
lap between the basis function fl^ (or ) around the 
point r n and the delta function around r k _i/ 2 and rfc+ 1 / 2 , 
which is the derivative of the basis function around r k , 


1 


drfln ( r ) 


dM k (r) 

dr 


/ 


drfln (r) 


dtt k '{r) 

dr 


= / drf2„(r)[d(r-r fe -i/ 2 )-<5(r-r fc+ i/ 2 )] =0, (15) 


drfln (r) 


,dV k (r) _ f dC]E > () dflf(r) 
dr / 71 U dr 


a (^fc,n+l ^k,n— l) 

2A 

(19) 

Then, the expansion of the fermion wave functions 

F(r)=^F„^(r), 

n 

(20) 

G(r) =^G n n»(r), 

n 

(21) 

uniquely discretizes the following terms in 
Hamiltonian 

the total 

J drF{r)G{r) 


= [ dr ^5jP»^(r)^fc(r)] 

J n k 


= EE^ G * / drfi'(r)^(r)] 

n k ** 



= ]T ^[F n G t Ad nfc ] = A ^2 F n G n , (22) 

n k n 


= I drfl%'(r)[5(r - r k _ 1/2 ) - S{r - r k+1/2 )] 


— $k,n— 1 H - (^fc,n $k,n) “i - ^fc,n+l 


for which we have used Eq. (HU- 

Similarly, using Eq. (fTiil) it follows that 


drF(r) 


dG{r) 

dr 


— ^fc,n+l $k,n— 1* 

(16) 

= / 


Using Eqs. flSJ) and ©, we further define 


J n k 


n nW = + O')]- 

(17) 

= EE[^g 1 a <4 ‘- +1 2 “ & --' ) ] 


With the help of Eqs. (ISll-dTTll. we get 




J drSl*(r)Sl%(r) 


A ITT Cr+1 ~ G ra _i 

2A 

n 

(23) 
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2.3 Matrix form of Dirac equation in 
terms of basis functions localized in 
spacetime 

In our scheme, the aforementioned total Hamiltonian in 
a fermion-confined region is expressed by 


n = 


A 


J2{[F n (+m- 


a 

r n 


crr n )F n 
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where for 1 < i, j < N, 

Q/ 

H Aij — Hij&ij ( -\-7Tl . ^ T O /A) 6j j , 


( 29 ) 


(30) 


r? + 1 G n - 1 TP ^ O' 

”n ~ A o n Or, 

2A r„ 


+ [G n {—m —-—h crr n )G n 

r n. 


o< F n+ 1 F n -1 o' K ip 

Gn 2A Gn v„ K 


- F n EF n - G n EG n }. 


Hx v = E v x v . 

A detailed form of the equation is 


Ha H b 
H c II n 


F" 

G 


= E v 


pv 

G' 


where the components of the row vectors x v , F 1 ’ and G 1 
have the relations 


x .V = R v 


V _ OIV 

x N+i 'Oi ) 


— K, 

iA 


HBij — — \ 


for j = i 
for j = i — 1 
for j = i + 1 
0 for the others 


+i 

2A 


-1 

2A 


(24) 


— K, 

iA 


H-Cij — HN+i,j — \ 


The potentials a/r n , or n and n/r n are replaced by 
a/(?iA), ntrA and K,/(nA), respectively. Let us impose 
the normalization condition on F and G and replace the 
last terms —F n EF n — G n EG n by —E(F n F n — 1/N) — 
E(G n G n — 1/AO, w hh AT being the number of lattice 
(grid) points, and E is regarded as a Lagrange multi¬ 
plier. 

By variational calculus, A remains in all the terms 
including those with E, and all terms are divided by A. 
We then get the Dirac equation in the matrix form for 
the w-th eigenvector x" and the associated eigenenergy 
E v 


for j = i 
for j = i — 1 
for j = i + 1 
0 for the others 


-l 

2A 


+ 1 
2A 


(31) 


(32) 


(25) 


(26) 


Hjoij — Hn — ( oi .^ T o/ A) Sjj. (33) 

3 Analytic eigenenergies of a 
fermion under a confining lin¬ 
ear potential 

3.1 Analytic eigenenergies of the Dirac 
equation with a linear potential 

The quantum eigenenergies are now treated by consid¬ 
ering the energies of a classical mechanical Hamiltonian 
for the principal properties of the Regge trajectory. We 
use the following matrix R in the same notation as the 
aforementioned matrix U 


R = 


(27) 

(28) 


where 


Rv 


Ra Rb 
R c Rd 


1 27 Tp r q r , 

a [sm( 2lvTT- 


(34) 


(35) 


for 1 < i < N. It is to be noted that i is not the imaginary 
unit in complex numbers, but an integer. The matrix H 
has the following form 


H = 


H a H b 
H c II i) 


for 1 < p r , q r < 2 N, and Gn is a normalization constant 
and large for large N. We then have (emphasizing the 
matrix elements of R using parentheses [ ]) 

(HR)i, qr = 

On 
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. 2ir{i)q r 

Fii[sm( 2 ivTT ) ] 


rr r • /27 t(TV + * — 1 )g r 

+ fi i ,jv +i -i[sm(- 2jv + 1 7 )] 


. rr r • , 2tt( 7V + *)<?,• 
+ [sm(- 


27V + 1 


, . /27r(iV + * + l)g r .. 

+ ff 1 , JV+ , + i[sm(-27VTl-( 36 ) 

The second and last terms below “1 /Cn{” in the above 
equation amount to an expression in term of the matrix 
component [sin(27r(TV + i)q r /(2N + 1))] of R 

lt r • r 27r ( N + i- l)g r 

^i,jv+i-i[sm(- 2N +1 -" 


, rr \ ■ ( 2 7r (-^ + * + 1 )Qr 


2iV+ 1 


1 . 2tt(N + i — 1 )g r 

2A Sm 2]V+1 ->1 


-1 . 2ir(N + i + 1 )q r 
2A Sm -2JVT1- 


-1 O 2^r(^ v+ »)gr \ 

J_ . / ~27T q r COS( 2JV+1 ) 

A sml 27V + 1 cm( 27r{N+i)qr ) 

2iV+l / 


. 2n(N + i)q r 
x [sm(- 


27V +1 ' J 
Thus, HR has been set to H'R with 


(37) 


H'i = H u , 


Hi-l,N+i — Oj 


(38) 

(39) 


Hi,N+i ~ Hi,N+i 


~ /2n( N-\-i)q 

1 . -27Tg r C08( \ N+ l q 

' A Sm 27V + 1 sini^zAAtili 

bm V 2AT+1 


■^r+l,AT+i 

where ft is contained in H l N+i . 
In a similar way, we obtain 

Hn+Lz-i = O ' 


(40) 

(41) 

(42) 


HN+i,i — HN+i,i 


1 . 2 nq r cos ( 2 ^{t) 

+ A Sm 47V+l j sin( ^0 r) > 

(43) 


(44) 

H'N+i,N+i = HN+i,N+i- 

(45) 


Therefore, the matrix H' has been reduced to 
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H’ = 
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X 

(46) 
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0 
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0 
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0 
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X 
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-.f 
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X 



Hamiltonian is Hermitian and R l = R, where R t is the 
transpose of R. The result for the above matrix H' es¬ 
sentially implies that HR = H’R = R t H t = R t H ,t . 
From H’ = R^H’ 1 ) and R t H n = HR, we get H’ = 
R{HR) = R t HR, which implies that this process is a 
unitary transformation. 

A set of four elements H' u , H[ N+1 , H' N+ii and 
H' N+i jv+i °f the above matrix is independent of the other 
matrix elements. The matrix composed of these four ma¬ 
trix elements is diagonalized using a unitary matrix 

{ u pq for p, q = i,N + i 

, (47) 

S pq for the others 

where u pq are the four elements of the (two-dimensional) 
unitary matrix. The determinant to yield the eigenener- 
gies is 

{E - H' U )(E - H' N+i N+i ) - H' i N+i H' N+i i = 0. (48) 

We now consider physically meaningful phenomena, 
whose eigenenergies are those for the row with q r = 1 
(the lowest oscillation case) of the aforementioned matrix 
R, and drop the ^-dependent term from the Hamiltonian 
owing to small | sin(±27rg r /(2TV-|-l))|. (The contribution 
of the string tension term aiA to the energy Hamilto¬ 
nian is large in the case of not small i within the strongly 
bound confinement regime, while n/(iA ) exhibits a sin¬ 
gularity for low i.) Additionally, we neglect the masses 
of the composite fundamental fermions, for the cases in 
which the masses are small compared to the rotational 
energy. We then get the eigenenergies 

A) 2 ^(M) 2 = 0, (49) 
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which states that 


a . ac 
E = -— + aiA+ ] —. 
i A iA 


(50) 


The above Coulomb term —a/(iA) is disregarded for the 
extremely small A (in some sense beyond the regime of 
usual computer simulations with the larger A) due to the 
asymptotic freedom of the non-Abelian field, yielding 


E = criA + 

lA 


(51) 


(i is not the complex number but an integer of the lattice 
index). 

Here, we consider a function of the number x in the 
real continuum in the region x > 0, defined by 

Ifcl Ifcl 

E x = —- + ax = (ax + -—). (52) 

x x 

The function E x above takes the minimum at 
x m = (M/cr) 1/2 as 


^min _ 


\k\ 


(MAx) 1/2 


f7(|K|/cr) 1/2 


= 2(1*1 a) 1/2 . (53) 

This x m is (considering N A > a: m ) measured with the 
lattice spacing A as 

x m = i m A + e E , (54) 

where i m is an integer and e E corresponds to a residual 
denoted as 

~ 7 - £E < T (55) 

In the limit as A —> 0, the residual e E vanishes (e E 0), 
and the eigenenergy E in Eq. (l5ll) at i = i m approaches 
the minimum value E m j n equal to E™ ln 7 giving 

Emin = (E™' n ) 2 = 4 |/ c |< 7 , ( 56 ) 

(which is independent of A). 

We note that |k| is the absolute value of the integer- 
type relativistic quantum number, which includes a sign 
and originates from the Dirac equation. It is noteworthy 
that only the secular equation gives rise to the absolute 
value. 

We compare our equality Eq. (15^1) given above with 
that obtained from another theoretical method [211 [23] 
using the classical mechanical Hamiltonian in spherical 
coordinates, which is briefly summarized as follows. The 
relativistic classical mechanical Hamiltonian comprising 
the kinetic energy and a linear potential is denoted as 

1 T (cl) = (P 2 + m 2 ) 1/2 + err, (57) 


where P is the relativistic momentum of a fundamen¬ 
tal particle. For small mass compared to the rotational 
energy, the classical mechanical Hamiltonian, which was 
described above, is reduced to 

i? (cl ) = P + ar. (58) 


Using the rotational quantity J, which roughly corre¬ 
sponds to |k| and which is written by 

J = Pr, (59) 

the classical mechanical Hamiltonian amounts to 

i7 (cl) = i +(T r. (60) 

r 

The energy minimum of H ( cl ) also occurs at 

r=( J ~) 1, \ (61) 

to give 

(^min) 2 = 4J(T. (62) 


The above relation essentially coincides with the afore¬ 
mentioned equality Eq. (1561) , which is consistent with the 
principal properties of the Regge trajectory [23J. 

From the experimentally observed slope 



0.93, 


(63) 


in natural units for the Regge trajectory [23J |28], the 
equality |k| = PA n /(4cr) results in l/(4cr) = 0.93 GeV -2 
to yield ^fo = 518.5 MeV. If we use the relation ^fa = 
2.255Amom, derived analytically in our previous paper 
[2j where Amom corresponds to the scale-invariant en¬ 
ergy of quantum chromodynamics (QCD), we arrive at 
an Amom of 229.9 MeV. This is larger than the 186 MeV 
calculated for a smaller setting of \fa = 420 MeV in our 
previous paper [2|. These values are consistent with the 
observed QCD scale-invariant energy of around 213 MeV 


3.2 Discussions 

Here, we add some discussions concerning the present ap¬ 
proach and results obtained in the previous section and 
subsection. Unlike our approach, the classical mechan¬ 
ical Hamiltonian with a linear potential and repulsive 
rotational potential [24l [25], described by Eqs. <57J- 
(|62D in this paper, has no basis in the Dirac/Yang-Mills 
equations. Furthermore, concerning the Regge trajec¬ 
tory expressed as a function of the rotational quantum 
number for eigenenergies (masses) of the pair of the con¬ 
stituent fundamental fermion and antifermion, the mech¬ 
anism yielding large binding energies compared with the 
composite fermion masses has not yet been fully clarified 
by other theoretical/numerical approaches. 

In contrast, we used the Dirac equation in spheri¬ 
cal coordinates at the first quantization level, considering 
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the following points. First, because spherical coordinates 
differ from orthogonal coordinates, the eigenenergies are 
derived as a function of the relativistic quantum number, 
considering the expression of the Regge trajectory. Sec¬ 
ond, the procedure for deriving the solution of the Dirac 
equation shows that the mass (binding eigenenergy) of 
the pair of the fundamental fermion and antifermion has 
its origin in the linear force between the particles and the 
angular potential proportional to the relativistic angular 
quantum number. Therefore, the eigenequation in the 
form of determinant gives rise to the lowest energy as a 
function of the relativistic quantum number correspond¬ 
ing to the classical mechanical Hamiltonian energy for 
the principal properties of the Regge trajectory [211 [25] , 
which is presented in the part containing Eqs. <53-© 
in this article. 

We note that our formalism for the Dirac equation 
in the present paper uses the linear potential with the 
attractive Coulomb potential, which was directly derived 
from the Yang-Mills equation using the path integral at 
the second quantized field-theoretic level. The poten¬ 
tial used here was calculated non-perturbatively using 
the Wilson loop, which has all orders of boson contri¬ 
butions. The Wilson loop Wq gave rise to (Coulomb 
potential + linear potential) in Eq. Q, as follows: 


V w (r) 


ln[W Q (r)] 
— t\ 


(64) 


where t 2 — t\ is the time interval used in calculating the 
Wilson loop. If we use the Wilson loop derived in our 
previous paper [2], the linear force part, which was pro¬ 
vided analytically, is expressed as 


Wql( r) = exp[-ar(t 2 - ii)], (65) 

where t\ is a small quantity. Then, the detailed form of 
Eq. when the rotational energy is large compared 

to the fermion mass, becomes 

E=V w (iA) + M (66) 

i A 

The eigenenergies of the Dirac equation obtained in this 
paper can be larger than the masses of the composite 
fermion and antifermion, when the masses are small com¬ 
pared to the rotational energy. The Polyakov line shows 
the deconfinement at high temperatures, and if we use 
the result of our previous paper [2], the Polyakov line P T 
which we analytically derived is expressed as 

P T = cos{arccos[exp(— arr)\ — arccos[exp(— ott £ )]} 


~ cos{arccos[exp(— arr)}}, (67) 

with r = 1 /(ksT) (fce and T are the Boltzmann con¬ 
stant and temperature, respectively) and r £ being a small 
quantity. Then, 

e , = - b,(P r ) « - ln[exp(=j2)] = A. = £1, (68) 


where e q is the binding energy of the pair of the funda¬ 
mental fermion and antifermion, and Eb = err. Equation 
(1681) shows that e q is small at high temperatures and the 
deconfinement of paired fermions occurs in some sense. 
The classical mechanical Hamiltonian describes the sys¬ 
tem at absolute zero (temperature) and did not treat this 
deconfinement at high temperatures. 

Before the next discussions, we briefly refer to the 
quenched case, which takes into account the Okubo- 
Zweig-Iizuka (OZI) rule [29-31], implying that the further 
fermion-antifermion pair creation is suppressed. (The 
quenched case corresponds to the following approxima¬ 
tion in the case of the path integral with respect to 
fermion Grassmann numbers T and T, including a ma¬ 
trix Mf, 


J dtE'd'F exp(—’FMf’F) = det(Mf), (69) 


which is set to unity.) 

In the operator formalism at Euclidian time t, the 
Green’s (two-point correlation) function is given by 

Gto(*)=<0|W B (t)^(0)|0>, (70) 

where < 0| and |0 > is the ground state vacuum and 
T~Lb is the Heisenberg-type Hamiltonian operator of the 
bound state, which is expressed in terms of the energy 
operator E as 

7is(i) = exp(Et)77B(0) exp(— Et). (71) 

Using Eqs. and <23, the s-th solution of the Dirac 
equation |s > with eigenenergy is denoted as 


|s >oc 



(72) 


(The above radial functions are multiplied by spin- 
angular components and a constant factor for the center 
of mass exp (iP^X^), with X^ and being the position 
coordinates for the center of mass and its momentum, re¬ 
spectively.) For the quenched case, the Green’s function 
given by Eq. <70j yields, using the above solution |s >, 


G t0 (t) =J2< °l#B(f)|a >< a|^(0)|0 > 

S 


= Y, I < °I^b(0)|s > I 2 exp(—£ (s) f), (73) 

S 

where the eigenenergy of the Dirac equation has appeared 
as a decay constant for the Euclidian time. 

Thus, the solutions of the Dirac equation associated 
with the eigenenergies enter into the operator formalism 
with the help of the OZI rule (the suppression of the 
further fermion-pair creation). As previously mentioned, 
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the Dirac equation indicates that the lowest mass (ex¬ 
pressed as an eigenenergy as a function of the relativistic 
quantum number when the rotational energy is larger 
than the constituent particle masses) of the pair of the 
fundamental fermion and antifermion originates in the 
confining linear potential and the angular potential pro¬ 
portional to the relativistic quantum number. The mass 
of the pair is allowed to be larger than the masses of the 
composite particles. Furthermore, the deconfinement fea¬ 
ture at high temperatures in some sense is described by 
the Polyakov line according to the Yang-Mills theory, as 
presented in our previous paper [2]. 

4 Conclusions 

We have presented a formalism for the Dirac field under a 
confining linear potential using basis functions localized 
in the spacetime continuum, which formulates fields of 
finite degrees of freedom. A given linear potential is that 
from the Wilson loop analysis for a non-Abelian Yang- 
Mills field. The Hamiltonian matrix has been analyti¬ 
cally diagonalized with the use of two sequential unitary 
transformations, thus yielding the eigenenergies of the 
confined fundamental fermion-antifermion pair. The low¬ 
est eigenenergy (as a function of the relativistic quantum 
number for the large rotation energy compared to the 
composite particle masses) is proportional to the string 
tension and the Dirac’s relativistic quantum number re¬ 
lated to the total angular momentum, which is consistent 
with the expectation. 
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